On computation of the total set of robust discrete-time PID controllers 
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Abstract — The problem of finding the set of all multi- 
model robust PID and three-term stabilizers for discrete-time 
systems is solved in this paper. The method uses the fact 
that decoupling of parameter space at singular frequencies is 
invariant under a linear transformation. The resulting stable 
regions are composed by convex polygonal slices. The design 
problem includes the assertion of intervals with stable polygons 
and the detection of stable polygons. This paper completes the 
solutions to both problems. 

I. INTRODUCTION 

It has been shown that the stabilizing region for time- 
continuous PID-controllers is defined by a set of convex 
polygonal slices normal to kp axis in the (kp, kj, kp>) 
parameter space. Different approaches which prove this 
result use generalization of the Hermite-Biehler theorem, see 
[2], calculation of the real-axis intersections of the Nyquist 
plot, see [3] and singular frequencies, see [6], [1] and [5]. 
In the cited articles, it has been shown that the design 
of PID parameters due to decoupling of PID controller 
parameters at singular frequencies may be divided into two 
subproblems: (A) assertion of stable intervals of parameter 
kp independently on parameters k\ and kp> and (B) detection 
of stable polygons on the plane (kj, kp>) for a given kp. the 
solutions to problems (A) and (B) for the continuous-time 
systems can be found in [1] and [7]. 

In the present paper the time-continuous theory is trans- 
ferred to the design of PID and three-term controllers in 
discrete -time domain. First related results are reported in [5], 
where the set of all Schur-stable controllers is shown to be 
composed by transformed polygonal slices in a parameter 
space (r\,r 2 ,r 3 ). While a solution to the problem (B) as 
proposed in [7] applies directly in the discrete-time domain, 



see Section IV this paper addresses the problem (A) in Sec- 



tion [V] Section III proves that parameter decoupling applies 
in general at singular frequencies. Basically, a function is 
constructed from the characteristic polynomial, such that its 
imaginary part on each point on unity circle depends only 
on the gridding parameter (here conventionally chosen) r 3 
and is independent on parameters r\ and r 2 . This result 
will lead to the conclusion that for a given plant a fixed 
number of singular frequencies must be available for its 
stabilization. Given that parameter r 3 defines uniquely the 
number of singular frequencies, one can directly discriminate 
the intervals of parameter r 3 , which for no triple (n, r 2 , r 3 ) 
can provide stabilization. The simplicity of the criterion is, 
however, achieved at the price of conservativeness, since 
just necessity is provided. To cope with conservativeness the 
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approach proposed in [1] for the continuous-time case may 
be directly used. All aspects of the method are illustrated by 
examples. 

II. Problem definition 

Consider a simple closed curve F = {z \ z = T(a) + 
j 77(a), a € [a, b]}, in z— plane, which is symmetric to 
the real axis z and can be expressed in the form 

T: F(t,tj)=0, (1) 

and a characteristic equation of the form 

p = A(z)Q(z,n,r 2 ,r 3 )+B(z)=0, (2) 

where A, B are polynomials in z and r\, r 2 , r 3 are real 
parameters, that enter linearly in Q, 



Q = 8 l (z)r l +8 2 (z)r 2 + 8 3 (z)r 3 , 



(3) 



where the polynomials 5, (z) are assumed to be of the order 
< 2 in z. 

The following problem is handled in this article: Compute 
the set of all parameters r\,r 2 ,r 3 , such that the polynomial 
Q is r~ stable, that is, all its eigenvalues are enclosed by F. 
Of main interest in this article are circles with center on the 
real axis T and an arbitrary radius, which will be referred to 
as T— circles. Especially important is the unity circle because 
of Schur-stability. T— region refers to the region enclosed by 

r. 

It may be easily shown that Q describes the characteristic 
equation of a feedback loop with a PID or a three-term 
controller. Indeed a discrete-time equivalent of the PID 
controller has the transfer function 

„ , ci + c 2 z + c 3 z 2 
(z+zi)(z-l) 

Its structure follows in the quasi-continuous consideration by 
applying the rectangular integration rule (s — > (z — \)/Tz) to 
the ideal PID controller kj/s + kp + kps, resulting in z\ = 0, 
or by the trapezoidal integration rule (s — > 2(z— l)/T(z+l)), 
resulting in z\ = 1. Also the realizable PID controller kj/s + 
kp + kjjs I '(1 + T\s) converts by the trapezoidal integration 
rule to the controller structure Q with a pole at z\ = 
— (27] —T)/(2T\ + T). The following derivation also holds 
for a three-term controller with an arbitrary second order 
denominator polynomial 

n(z){c\ +c 2 z + c 3 z 2 ) 



Ci{z) 



d{z) 



(5) 



For both controller structures, Q and Q, the polynomial Q 
is of the form 



Q = c\ +c 2 z + c 3 z 



(6) 



It is clear that ([3]) and ([6]) are connected via a parameter 
space transformation 



c = T r, det T ^ 0, 



(7) 



with c = [ci,C2,ct,] t , r = [r\,r 2 ,r 3 ] . The rotation matrix T 
is determined by the polynomials, 8i(z),5 2 (z),8$(z). 

III. Basic definitions and theorems 

Let H and G represent the real and imaginary part of the 
characteristic polynomial p{z,r\,r 2 ,r 3 ) in ^ on T. 

Definition 1: F is said to be singular with respect to 
parameters r\ and r 2 in ([2| 



Rank er) = 1 



for all z G T. 



(8) 



d{r\,r 2 ) 

The latter equation is referred to as the rank-condition. 
Unless otherwise stated, F will be assumed singular for the 
rest of the article. A zero of the polynomial which addi- 
tionally fulfills the rank-condition is referred to as singular 
frequency. For a fixed parameter r 3 , singular frequencies 
are mapped to straight lines on the plane (r\,r 2 ) and the 
eigenvalues of Q can cross over a singular F just at singular 
frequencies. 

Definition 2: A function Ey{z) defined as 



Q(z,n,r 2 ,r3)=Er(z) q(z,r h r 2 ,n) 



with 



I q = ng-i{a) +g (a), a e [a, b] 



(9) 



(10) 



i.e. — = 0, vv/iere stands for the imaginary part of 
q, will be referred to as decoupling function of Q over F. 
Lemma 1: IfY is singular then 

dla dl a 

^=o^^=o, Vzer. (li) 

dr\ dr 2 
Proof. Consider the rank-condition d8). It may be checked 

that 



d(H,G) 
d(r u r 2 ) 



{R 2 a+4)(R 



2 



-II) 



dR q dl q 



dl q dR q 



dr\ dr 2 dr\ dr 2 



If rank-condition applies, that is, the rig^ht-hand side of the 
above equation is zero, then = =>■ -gf- = 0, since ^ 

0. Similarly, |^=0^^=0. □ 
Theorem 1: Two decoupling functions of Q in (JiJ over T 
are simply given by 



E T {z) = 8 x {z), £r(z) = 52(z) 
Proof. Suppose £r(z) = <h(z). Then 



(12) 



q = n 



Si (z) , fis(z) 



-'"2 



5i(z) 



-'"3- 



a/ 



and t- 1 = 0. Lema 



guarantees that for all z G T the 
imaginary part of is also zero, that is, ( jioj i applies. 



Since 



52 U) 
5i(z) 



is real on T, its inverse will be also real, that is, 



the imaginary part of |4§j for all z G T is zero. Hence, ^(z) 
represents also a decoupling function. □ 
The next statement may be directly checked. 



Theorem 2: Consider the function 

P{z) 



F(z): = 



A(z)E r (z) 



(13) 



The equation F(z) =0 for z G T decouples the parameters 
r\,r 2 and r$ into two equations 

rihi(a)+r 2 h 2 {a)+h Q {a) = 0, (14) 

r 3gl {a)+g {a) = 0. (15) 
The latter theorem is essential for solving the problem (A) 

as defined in the introduction of the article. The first equation 

refers to the generator of singular lines, and the second one 

to the generator of singular frequencies. 

Finally, let {p v } represent a finite set of polynomials of 

the form |2]), i.e. 



Pv =A v (z)Q(z,n,r 2 ,r 3 )+B v (z). 



(16) 



For example, this may be a multi-model of a continuum of 
plants or Kharitonov polynomials of an interval uncertainty. 

Theorem 3: The rank-condition for the set of polyno- 



mials (16) does not depend on the polynomials A v (z) and 

B V (Z). 

Hence a singular F is completely determined by the 
polynomial Q in Though the singularity of F is not 
destroyed by polynomials A v (z) and B v (z), they influence 
the root-condition, i.e. the singular frequencies on F. 

A. Schur- stability 

Consider the Schur-circle defined as F[ = {e^ a : a G 
[—71, 7C]}. It can be easily checked that the rank-condition 
|8j is satisfied over F{ for 



Q= (l+z 2 )ri+zr 2 + r 3 
and that the matrix T, as defined in |7| is 

10 1 

1 

1 



(17) 



(18) 



According to Theorem [T] the decoupling functions of ( 17 1 
over the Schur-circle F\ are 



E r (z) = l+z 2 or E r {z) 



(19) 



Let Ey(z) = z. It may be easily shown that the imaginary 
part of the function 

F(z) = —r l +r 2 + B j^r 3 (20) 
z zA(z) 

on the Schur-circle is of the form ( fl~0] >, that is imaginary part 



of 



l+z 4 



is zero on the Schur-circle. 



B. F-stability 

Consider a T— circle with center on real axis F = {m + 
re ;a , a G [—It, 7t]}, Fig[T] Now it can be shown that 



Q= {r 2 -m 2 +z 2 )r\ +(z-m)r 2 + r 3 . 



(21) 



For circles with center at x =m and radius r a transfer- where C a = cos(a). For illustration purposes its plot is 
mation matrix from c— to r— parameter space is checked to shown in Fig. [2] E.g. for r3 = —0.26118, the singular 



be. 



T = 



r — m 

1 



—m 1 
1 




A corresponding decoupling function is 

Er(z) =z — m 



and 



F{z) 



-m 2 +z 2 



z — m 



-r\ + r 2 - 



B(z) 



(z — m)A(z) 



l~3- 



(22) 



(23) 



(24) 




frequencies lying on the Schur-circle are computed to be 



a[=0 

a' 2 = ±0.4097 
a' 3 = ±0.97300 
a' 4 = ±n 



zi = l, 

z' 2 =0.9172 ± jO.3983, 
4 =0.5628 ± ;0.8266, 
z' 4 = -l. 



(30) 




0.5 



Fig. 2. The generator fo singular frequencies; Equation |29| 



Fig. 1. Schur- and I — circles 

C. Example 

Consider the discrete-time model of the plant 



G= 10 



_ 6 4. 1 65z 3 + 45 .77z 2 + 45 .77z + 4. 1 65 



z 4 -3.985z 3 +5.97z 2 -3.985z+l ' 
and a three-term stabilizer 

- 2 - 1.54U + 0.5992) (ei + c 2 z + c 3 z 2 



C(z) = 10' 



z(z + 0.4047) (z + 0.2162)(z - 0.4934) 



(25) 



(26) 



whose parameters c\,C2,c 3 should be synthesized. The syn- 
thesis is done in (ri,r2, ^-parameter space. Therefore the 
transformation (flBl can be used. Then 



A = z 5 +9.44z 4 -5.34z 3 -9.34z 2 ±5.04z + 0.59(27) 
B = 0.19z 8 -0.73z 7 +z 6 -0.45z 5 -0.12z 4 + --- 



0.14z -0.009z 2 -0.008z. 



(28) 



The equation ( 15 i reads 

C 5 a - 2.7 1 Ci + 1 .61 Cl + 1 .35 C 2 - 1 .70C a + 0.44 



0.02C 4 + 0.08 Cl - 0. 1 3 C 2 - 0.08 C a + 0. 1 1 



(29) 



IV. Stable convex polygons 

Notice that the characteristic equation in Q over singular 
r decouples equivalently to the two equations in ( 14 1 and 



(15 1. For a fixed rj the equation ( 14 1 defines a set of singular 



frequencies {cc' v } on T. The equation (15i determines the 
set of straight lines {A v }, which are determined by {«(,}. 
Geometrically, the straight lines {A v } represent the eigen- 
value boundaries in the (ri,r2)-plane for the fixed parameter 
r3. Thus, stable regions are composed by convex polygons 
and the design problem defined in the previous section is 
decoupled into two subproblems: (A) assertion of stable 
intervals of parameter r3 independently on parameters r\ and 
r 2 and (B) detection of stable polygons on the plane {r\ , r 2 ) 
for a given r^. 

For the sake of completion, this section recalls briefly 
the solution of problem (B), while problem (A) will be 
thoroughly discussed in the next section. For details the 
reader is referred to [7], [1]. The algorithm is motivated 
by the concept of inner polygons, which claim a necessary 
condition for stability: A polygon n is said to be an inner 
polygon if any transition over X v inside the polygon causes 
an eigenvalue-pair to enter the T— region at the corresponding 
frequency. 

In order to automate the detection of an inner polygon 
each singular line X (a 1 ), will be assigned a "transition" 




To 




Tl 



Fig. 3. Definition of e\ and ey. the motion of eigenvalues in the vicinity of 
a singular frequency z'. By convention, the shaded part refers to the stable 
side of r and the normal N points outside the T-region. 



function e: it is negative if the transition [8r\ , 8r2\ (see Fig. [3]) 
over the singular line causes an eigenvalue to become stable, 
otherwise it is positive. Let e\ correspond to 8r\ > 0, 8r2 = 
and <?2 to 8r2 > 0, 8r\ = 0. In [7] it shown that the motion 
of eigenvalues at z' = T(of')± jri(a') defined as the scalar 
product e = jj,N (see Fig. [3]l is described by the equation 



«i/2 - 



/ dF 


d(H,G) 


dF 


d(H,G) 






drj 





(31) 



Using this information, an algorithm for the detection of the 
inner polygons with maximal T-stable eigenvalues may be 
developed, see [7]. Such polygons are the only candidates, 
which need to be checked for stability. To this end, it suffices 
to check any point within such a polygon. 




1 .8 2 2.2 



Fig. 4. The stable polygon lying on the plane = —0.26118 

A. Example: (cont.) 

Reconsider the plane r?, = —0.26118. The corresponding 
singular frequencies were computed in the previous section, 
(30 1. Each of these singular frequencies z' v generates a 
straight line X v on the plane tt, = —0.26118. The resulting 
stable polygon is shown in Fig. [4] It is enclosed by the 
straight lines frequencies X\,X% and A3, corresponding to 
z',,4 and z' 3 . 



V. Stable gridding intervals 

This section focuses on the problem how to discriminate 
T3— intervals, such that stable polygons may exist therein. 
The basic idea presented here tempts to extrapolate this 
information from the plot of equation (15 1, see Fig. [2] 
Basically, one would like to link somehow the stability of 
the characteristic polynomial Q with the number of singular 
frequencies on Y. This approach is motivated by the fact 
that an emerging pair of singular frequencies produces new 
polygons and vice-versa. 

The following lemma is essential for the main result 
of this section. Without loss of generality, the T-region 
is assumed to be a Schur-circle. The generalizations for 
T— circles are straightforward. Its proof is, however, beyond 
teh scope of this article (proof hints: (a) N is always even, 

(b) the Mikhailov curve is symmetric to the real axis T and 

(c) the poles on T are avoided by infinitely small circles, 
whereby the Mikhailov plot experiences a phase shift of ±7T 
at infinity.) 

Lemma 2: Consider the Mikhailov plot of a function 
F(e ja ) on the Schur-circle T\ and let F(z) possess J poles 
=/= ±1 on T\, a pole of the order J + at z — +1, and a pole of 
the order J at z = — I. If the phase change of the Mikhailov 
vector F(e jCC ) over T — 1 is N% as a changes from to +7Z 
without touching them, then it cuts the real axis at least 
Z— times, where 

N-J-2-E(J+)-E(J) 
Z= — (32) 

Theorem 4: Consider the characteristic polynomial (|2]) 
and the Schur-circle T\, Let 
N: order of the polynomial A2u 
R: number of zeros of A(z)Er(z) lying inside T 
J: number of zeros 7^ ±1 of A(z)Ey{z) lying on T 
J + : order of the zero +1 of A(z)Ey{z) 
J : order of the zero —1 of A (z)£r(z) 

Z: number of singular frequencies in the interval < a < 

is 



A necessary condition for stability of |2]) 

J + E(J + )+E(J) +2 

Z>N-R — . 

2 

Proof. As shown in Section [Hi] the equation 

P(z) 



F(z) = 



= 



(33) 



(34) 



A(z)E r (z) 

decouples parameters r\^r% and rj into the two equations 



(14 1 and (15 1, whereby the imaginary part represents the 
generator of singular frequencies. The Mikhailov plot of 
function F(z), z = t(a) + jrj(a)), - % < a < % intersects 
the real axis exactly at singular frequencies z' =z{ot'). If p(z) 
is T— stable, then according to the principle of argument on 
T\ yields 

A4> F = {N-R)2n, (35) 

where A<Pf represents the phase change of the function F(z) 
on F\ as a changes from — n to +7T, including ±n. Notice 
that J + J+ +J- zeros of A(z)Er(z) lying on T are avoided 
by exclusion via infinitely small semicircles, i.e. they are 



considered to be outside Tj . Thus according to Lema [2 the 
equation (33 i results if N is substituted by 2(N — R) in (32i. 

□ 

Using this theorem one can directly read from the plot 
of equation (15i the — interval(s) where stable polygons 



may exist. However, the bounds defined by this theorem 
may be conservative if stable polygons disappear due to the 
intersection of at least three singular lines at one point in 
the parameter space (n,^,^). The reader is referred to [1], 
where a thorough discussion on this item is provided. 

A. Example: (cont) 



Consider Schur-stability for A and B given in ( 27 1 and 



(28 i. It can be checked that A(z) possesses three zeros inside 
the Schur-circle, one zero at z = — 1 and one zero outside the 
Schur-circle. Thus if the decoupling function £r(z) = Z, see 
d23j, is used, it follows that 

N = 8, # = 3 + 1=4, 7=1, /+=0, and J = 1. 

Hence, for stability, Z > 3 singular frequencies are required 
in the interval < a < +n. In order to discriminate stable 
rj intervals one should check the plot of the generator of 
singular frequencies shown in Fig. [2] The stable interval is 
depicted by the grayed strip in Fig. [2] 

-0.52236 < r 3 < 0.00290. 

Notice that the zoomed plot in Fig. [2] recommends that for 
< r3 < 0.00290, four additional singular frequencies appear. 

On the other side if the decoupling function Er(z) = 1 +z 2 
is used, then 



N=l 



R = 3, 7 = 3, /+ = 0, and J = 1 



i.e. again for stability Z > 3 singular frequencies are required 
in the interval < a < +n. 

B. Example: PID control 

Now consider PID control of the same plant ( |25] > with the 
control law Q. It can be shown that in this case, 

A = z 3 + 10.98z 2 + 10.98z+l 



B = 0.1z 6 -0.5z 5 



(36) 

-0.5z 2 -0.1z. (37) 

By using Ey(z) = z, it is easily checked that 

N = 6, R= 1 + 1=2, / = 0, /+ = 0, and /_ = 1 

since A(z) has one zero inside the Schur-circle, one zero at 
z = — 1, and the third one outside. Hence, Z > 3 singular 
frequencies within < a < +n are required. However, one 
can check that the maximal number of singular frequencies 
within < a < +n is 2, so no PID controller can stabilize 
the plant ((25]). 

C. Example: (cont) 

Having identified the stable intervals, griding of parameter 
r$ yields stable slices, which compose the 3-D stability 
region. E.g. Fig. [5] shows the whole region in (n,^,^) — 
parameter space of three-term controllers, (26 1, which stabi- 
lize the plant p5). 




Fig. 5. The set of all Schur-stable polygons 



VI. Conclusion 

The problem of finding the set of all PID and three-term 
stabilizers for a linear discrete-time system is treated in this 
paper. The basic result is the transfer and generalization of 
the counterpart theory of continuous -time PID stabilizers to 
time-discrete domain. The design method is based on the fact 
that controller parameters in a linearly transformed parameter 
space appear decoupled at singular frequencies. Thereby the 
design problem decouples into detection of stable polygons 
and assertion of intervals where such polygons may exist. 
A new simple and powerful rule is introduced to discrim- 
inate such intervals. The design method for simultaneous 
stabilization of several operating points becomes feasible by 
intersecting convex polygons. 
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